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Abstract 

Dynamical systems methods are used to investigate domain-wall solutions 
of a two-parameter family of models in which gravity is coupled to an axion, 
and to a dilaton with an exponential potential of either sign. A complete global 
analysis is presented for (i) constant axion and (ii) flat walls, including a study 
of bifurcations and a new exact domain-wall solution with non-constant axion. 
We reconsider 'fake supergravity' issues in light of these results. We show, 
by example, how domain walls determine multi-valued superpotentials that 
branch at stationary points that are not stationary points of the potential, and 
we apply this result to potentials with anti-de Sitter vacua. We also show by 
example that 'adapted' truncation to a single-scalar model may be inconsistent, 
and we propose a 'generalized' fake supergravity formalism that applies in some 
such cases. 
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1 Introduction 

In a recent paper pQ we addressed the problem of domain wall solutions of the cou- 
pled Einstein-dilaton equations, in d spacetime dimensions, using dynamical systems 
methods imported from studies of cosmological solutions of the same model [21 E]. 
We recovered efficiently and simply many of the previously known exact results on 
dilaton domain walls and found some new ones. In addition, we obtained a qualita- 
tive overview of the entire space of domain wall solutions in this model. For example, 
the "Janus" walls that were studied in |5] were shown to have a natural interpreta- 
tion as marginal bound states of a new type of dilaton domain wall that we called 
a "separatrix-wall" . This result is reminiscent of BPS solitons and it suggested a 
'hidden' supersymmetry of the separatrix-wall that has since been confirmed [5]. 
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One purpose of this paper is to continue our dynamical systems investigations of 
domain walls in a more general class of models for which both the dilaton field a 
and a pseudo-scalar 'axion' field x are coupled to a <i-dimensional metric g via the 
Lagrangian density 



where \i is an axion-dilaton coupling constant, and the potential function V takes the 
form 



where A is the dilaton 'coupling constant', and A a non-zero 'cosmological' constant. 
This is the unique form of the potential that preserves an invariance of the equations 
of motion under the 'dilation' g — > e w g, for constant parameter u, by virtue of the 
inhomogeneous transformation a + to/ A of the dilaton field. Only the sign of A is 
relevant when A 7^ 0, as a shift of a is then equivalent to a scaling of A. We may 
assume that A > without loss of generality, as in p], but then \i could have either 
sign; its absolute value is related to the radius of curvature of the hyperbolic target 
space for which (a, x) are coordinates. The \i = case, which corresponds to a flat 
target space, has been investigated (in a cosmology context) in J6j[7]. 

We have recently studied cosmological solutions for this class of axion-dilaton 
model [8], following the work of [9] on a particular model for which flat universes 
were shown to expand and contract quasi-periodically in a certain non-Einstein frame. 
We found a range of the parameters (A,/i), for which generic flat universes are eter- 
nally expanding, in Einstein frame, but undergo a medium-time and/or late-time 
oscillation between acceleration and deceleration. These results were one motivation 
for the present work because if one allows for both signs of the potential then the 
family of dynamical systems governing cosmological solutions is the same as the fam- 
ily of dynamical systems governing domain-wall solutions (illustrating the general 
'Domain- Wall/Cosmo logy Correspondence' [5]), and the quasi-cyclic nature of many 
trajectories could imply novel behaviour for associated holographic renormalization 
group flows. 

The equations for either domain walls or (homogeneous and isotropic) cosmologies 
in a model with n scalar fields define an autonomous dynamical system of dimension 
2n + 1, once account is taken of the reparametrization constraint. Reparametrization 
invariance implies that the physical phase space is actually only 2n-dimensional but 
the dynamics on this space is not necessarily autonomous. Thus, even for n = 1 
we should expect to have to consider a 3-dimensional system, and a 5-dimensional 
system for n = 2. However, a reduction to an autonomous 2-dimensional system is 
possible for n = 1 when the potential is an exponential, as is well-known; this case 
arises in the present context from the restriction to constant axion. Less well-known 
is the fact that the 5-dimensional dynamical system for n = 2 models of the above 
axion-dilaton type has a 2-dimensional autonomous subsystem when the restriction is 
made to flat domain walls, or cosmologies, and the behaviour is essentially determined 
by this subsystem. 




(d X f ~ V 




V = Ae 
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Here we present a detailed study of the 2-dimensional dynamical systems corre- 
sponding to (i) constant axion and (ii) flat domain walls, for either sign of A and 
as a function of the two parameters (A,//). Dynamical systems methods are partic- 
ularly powerful in these cases because trajectories are easily visualized, and chaotic 
dynamics is excluded. Case (i) has been investigated previously in [2], [3], [7] (in the 
context of cosmology) and in our previous paper [1] (in the context of domain walls) 
but here we present a complete global analysis. In particular, we show that the global 
phase space, allowing for both signs of A, is a sphere. We also extend the analysis 
of bifurcations to those that occur 'at infinity', which corresponds to a great circle 
on the sphere separating the trajectories with A > from those with A < 0. Case 
(ii) has been investigated previously in the context of cosmology [U [9] for one sign 
of A. Here we present a complete global analysis that includes both signs of A. One 
new result is an exact flat domain- wall 'scaling' solution with non-constant axion 
field. It corresponds to a fixed point of the dynamical system and for A < it is the 
domain- wall analog of the A > axion-dilaton cosmological solution found in [8]. 

For some values of the parameters (A, /i) and a choice of the sign of A, the model 
defined by (jl.ip may be a truncation of a supergravity theory. For example, the 
Freedman-Schwarz N — 4, d — 4, supergravity [TU] has 2A = — fi = 2, ancjj A < 0, 
and in this context one may ask whether any given domain wall solution preserves 
some fraction of the supersymmetry. A necessary condition for partial preservation 
of supersymmetry is that the domain wall solution admit a Killing spinor, which is 
a spinor field satisfying the equation that results from the requirement of vanishing 
supersymmetry variation of the gravitino field(s); this equation depends not only on 
the spacetime metric but also on the (pseudo) scalar fields through a 'superpotential', 
which determines the scalar potential. More generally, one can extend the notion of a 
'supersymmetric' domain wall solution beyond the supergravity context by allowing 
the superpotential to be any function that yields the potential according to a d- 
dimensional generalization of the formula that applies in the supergravity case 
Domain wall solutions admitting Killing spinors with respect to such a superpotential 
[T2J, [13] are then said to be 'fake supersymmetric' solutions of a 'fake supergravity' [I]. 
The relevance of this idea is that 'fake' supersymmetry suffices for classical stability. 

Another purpose of this paper is to use the axion-dilaton models as a 'laboratory' 
for further investigations of fake supersymmetry, although we will also have more 
to say about the simpler Einstein-dilaton model of pQ. It was shown in [U [T] that 
'almost all' flat domain-wall solutions of any single-scalar model are fake supersym- 
metric because the solution can be used to construct a superpotential with respect to 
which the first-order Killing-spinor integrability equations are satisfied. This result 
was extended to particular curved domain walls in [13] and to 'almost' all curved 
domain walls in [5]. The qualification 'almost' arises from a monotonicity condition. 
Typically, this condition is violated at isolated points (although non-isolated accu- 
mulation points occur in walls that are asymptotic to perturbatively unstable anti-de 
Sitter vacua [H]). In this case the associated superpotential is multi- valued and the 
domain wall is 'piecewise' supersymmetric with respect to it [5]. We illustrate this 

1 This model was incorrectly identified in [5] as one with A > 0. 
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phenomenon here with a simple double- valued superpotential that can be found from 
an exact flat domain wall solution of pQ. This example illustrates another significant 
feature, intrinsic to the phenomenon: the branch point of the superpotential is a 
stationary point that is not a stationary point of the potential. We use the insights 
gained from this example to discuss general potentials with adS vacua, showing that 
there again exist conditions under which some superpotential will be double-valued 
and hence not defined for all values of the scalar field, as was recently argued for a 
particular potential [15] . 

Fake supergravity models with multiple scalar fields were investigated briefly in 
[T3] and more extensively in [16j, where it was observed that for any given domain- wall 
solution there exist, at least locally, 'adapted' target space coordinates such that all 
scalar fields but one are constant [IE] . The given solution is thus manifestly a solution 
of the single-scalar model obtained by the obvious truncation in the new variables. 
It is also (if it satisfies the monotonicity condition) a supersymmetric solution of the 
'adapted' single-scalar model [5]. Thus, results on fake supersymmetry for single- 
scalar models can be extended simply to multi-scalar models. However, there is a 
difficulty with this extension (in addition to possible global problems arising from the 
local nature of the 'adapted' target space coordinates). To address stability of a given 
solution one needs to consider solutions 'nearby' the given one. It may happen, if 
the single-scalar truncation is inconsistent, that there are no nearby solutions of the 
single-scalar model that are also solutions of the multi-scalar model, in which case 
the fake supersymmetry of the domain wall as a solution of the single-scalar model is 
irrelevant to its stability. The consistency of the single-scalar truncation is therefore 
important. Consistency is manifest for the restricted class of multi-scalar models 
considered in [13] and it was shown in [16] that one-scalar truncations adapted to 
supersymmetric do main- wall solutions of d = 5 supergravity are consistent. 

In view of these facts, it would be a rather unsatisfactory state of affairs, violat- 
ing the sprit of fake supergravity, if a truncation adapted to a non- supersymmetric 
solution of some supergravity theory were to be consistent, so it is natural to con- 
jecture that this never happens, i.e. that any one-scalar truncation adapted to any 
non- super symmetric domain-wall solution of some supergravity theory is inconsistent. 
For a solution that is not known explicitly we know of no way to determine whether 
the associated adapted truncation is consistent, so we are not in a position to prove 
the conjecture. However, the new exact domain-wall solution with non-constant ax- 
ion mentioned earlier allows us to test it. The target space coordinates adapted to 
this solution can be found explicitly and it turns out that the associated single-scalar 
truncation is inconsistent. If the range of the parameters (A, fi) for which this solution 
exists are compatible with the consistent embedding of the model in a supergravity 
theory then the conjecture states that the solution must be non-supersymmetric in 
this context. As we shall show, there is a consistent truncation of minimal d = 4 
supergravity to the d = 4 axion-dilaton model provided that A/x = —2, and this re- 
striction on the parameters is compatible with the existence of the new exact domain- 
wall solution, which may therefore be viewed as a supergravity solution. We then 
show, in agreement with the above conjecture, that it is not supersymmetric in this 
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context. Moreover, it remains non-supersymmetric as a solution of a <i-dimensional 
axion-dilaton fake supergravity that we propose as a generalization of the standard 
(i-dimensional one-scalar fake supergravity formalism. 



1.1 Preliminaries 

As in p, we introduce the (^-dependent constants 



d- 1 



" = 1/2^, j3 = l/y/2{d-l){d-2). (1.3) 

and express the spacetime metric in terms of a function (p(z) through the domain- wall 
ansatz 

ds 2 = e 2 ^f 2 dz 2 + e 20v dZ* , (1.4) 

where f(z) is an additional arbitrary function that allows us to maintain invariance 
under z-reparametrizations, and dEf, is the metric of a (d — l)-dimensional maximally 
symmetric spacetime with inverse radius of curvature k; i.e. such that the scalar 
curvature is k(d—l)(d—2). We may restrict k to take the values —1, 0, 1, and we then 
have a domain wall spacetime that is foliated by anti-de Sitter (adS), Minkowski or de 
Sitter (dS) spacetimes, respectively. Invariance of the solution under the isometries 
of the metric implies that the fields (a, x) can depend only on z. The field equations 
now reduce to equations of motion, and a constraint, for the variables (<£>, a, x) that 
are the Euler-Lagrange equations of the effective Lagrangian^] 

Left = \ if' 1 (V 2 -° 2 - e^X 2 ) + /e 2 ^ {k(3- 2 e~ 2 ^ - 2A e ^)] . (1.5) 
The equations of motion for the gauge choice 

are equivalent to the equations 

Cp = -a^+ 1 ~\pa-2aA+J^e^, 

a = -aa^p + - (A - p) a 2 + i/i0 2 + (// - A) A (1.7) 



and the constraint 



X 2 = 



p 



'U 



The x equation of motion follows from differentiation of the constraint. 



2 Because of the inclusion of the function / in the ansatz, this effective Lagrangian can also be 
obtained by substitution of the ansatz in the Lagrangian density (jl.ip , whereas this is true only for 
k = for an ansatz that fixes the z-reparametrization invariance. 



5 



Although the gauge choice (11. 6p is convenient for our particular models, it has the 
disadvantage that z is not (generically) an affine parameter when A 7^ 0. It is related 
to an affine distance parameter z through the differential equation 

dz = e^ Xa{z) dz, (1.9) 

and the spacetime metric takes a standard form in terms of z; for flat walls we have 

ds 2 = d~z 2 + e 2 ^ds 2 (E^ d -V) , (1.10) 

where (perpetrating a slight abuse of notation to avoid the introduction of new sym- 
bols for dependent variables) <p is now to be understood as that function of z obtained 
by use of the inverse function z(z) in <f(z), and similarly for a. 

For the two special cases (i) % = and (ii) k — 0, the equations (11.71) contain a 
2-dimensional autonomous dynamical system for the variables 

u = & , v = ip . (1-11) 

In the first case, two values of A, the 'critical' value A c and the 'hypercritical' value 
A/j are special because they correspond to bifurcation points. In our conventions 



\ c = 2^= ^J — -, \ h = 2a = Vd^T\ c . (1.12) 

2 Constant axion 

For x — 0, the equations (11.71) and constraint ( 11. 8ft yield the autonomous dynamical 
system 

ii = —auv + 2^ u2 ~ ^ 

v = —(5v 2 u 2 H — Aw — 2/3A (2.1) 

and the constraint 

v*-u 2 + 2A=-^e*- 2 to. (2.2) 

P 

The constraint shows that the k — — 1 trajectories are separated from the k = 1 
trajectories by the k = hyperbola v 2 — u 2 + 2A = 0. As the absolute value of A 
is irrelevant for A > 0, we effectively have two one-parameter families of dynamical 
systems, one for A > and another for A < 0. In the absence of limit cycles, 
the phase-plane portraits of these two systems are determined by the fixed points. 
Poincare index considerations show that limit cycles cannot arise for A > 0, nor for 
k < [lj. This observation leaves open the possibility of k = 1 limit cycles in the 
A < system, but in this case the constraint ( 12.21) implies that v 2 > 0, and hence 
that tp{z) is a 'monotone function', which is incompatible with the existence of a limit 
cycle. 
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We now present an analysis of the global structure of the A > and A < 
dynamical systems. As we shall show, individually the global phase space for each 
system is a disc, but the two discs can be viewed as two hemispheres of a topological 
spheral 



2.1 Global structure 

Define new phase-plane coordinates (x+, y + ) by 



{x+,y+) 



l 



Vu 2 + v 2 + 2A 



[u,v) 



(u,v) 



2A 



2 y~\ 



\-x\-y 



(2.3) 



This maps the k = hyperbola u 2 — v 2 = 2 A into the lines 2x\ = 1. For A > the 
entire (u,v) plane is mapped to the interior of the unit circle in the (x+,y+) plane. 
For A < the domain of the map is the exterior of the circle of radius a/2|A| in the 
(u, v) plane, which is mapped to the exterior of the unit circle in the (x, y) plane. 
Thus, the new coordinates yield a description of the dynamical system that includes 
both signs of A, with A effectively changing sign on the unit circle; this is possible 
because the value of A becomes irrelevant as the unit circle is approached. 
The equations for the dynamical system in the (x+,y+) coordinates are 



dx + 
dz + 
dy+ 



1 



x+y+ [x\ - yl) + (x 2 + y 2 - l) 



~ 3y _ 



(4 - yl) + (4 + y 2 + -l) 



-A — AX i + — S + Uj 

(3 - \x + y + + ^-x 2 , 
Aft 



(2.4) 



where z + is a new independent variable such that 



2A 



l-x%-yl 



(2.5) 



Observe that the unit circle in the (x + , y + ) plane is a fixed set of the system defined by 
the equations (12.41) . Defining polar coordinates (r, 9) by x + = r cos 9 and y + = r sin 9, 
we find that 

r'l =0, 9'\ = -— cos#cos2#. (2.6) 

r=l r=l \ v 7 

Aft 

This determines the flows on the unit circle, on which there are six fixed points, 
two at the angles for which cos# = and four at the angles for which cos 2(9 = 0. 
The trajectories for r < 1 are those of the original dynamical system restricted to 
A > 0. Similarly, the trajectories for r > 1 are those of the original dynamical system 
restricted to A < 0, but the description of this system in the (x+,y+) coordinates is 



3 This result was stated without proof in [S], where phase planes for the two hemispheres were 
sketched for some cases of relevance to cosmic acceleration. Only A > cosmologies were considered 
in [8], since cosmic acceleration is possible only in this case, but it should be remembered that this 
corresponds to A < in the domain wall context, and that the sign of k is different. 
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incomplete because the region in the (u,v) plane with u 2 + v 2 < 2|A| is not covered 
when A < 0. 

We may remedy the incompleteness of the (x + ,y + ) coordinates for A < by 
defining a different set of new coordinates by 



(x-,y-) 



1 



Vu 2 + v 2 - 2A 



(u,v), (u,v) 



2A 



x 2 ^ + yi - 1 



(*_,!/_). (2.7) 



This maps the fc = hyperbola u 2 



2A into the lines 2y 2 _ = 1. The domain of 



this map for A < is the entire (u, v) plane, which is now mapped to the interior of 
the unit circle in the (x_, y_) plane. It is now for A > that the domain is the region 
in the (u, v) plane exterior to the circle of radius v^2A, and this region is mapped to 
the exterior of the unit circle in the (x_,y_) plane. The equations for the dynamical 
system in the coordinates are 



d/X — X / 2 2\ / 2 2 i\ 

j — = T-aJ-J/- - Z/-J + + V- ~ 1J 

a.z_ Aft 

^ = -f* 2 _(* 2 _-y!) + (*_2 + y!-l) 
az_ Aft 

where Z- is a new independent variable such that 



-/3 + /3y 2 - + 



Aft 



(2.8) 



2A 



2;_ = 



(2.9) 



The unit circle is again an invariant set of the system, and the dynamics on this set 
is exactly the same as it was for the (x + ,y + ) coordinates. 

We now have two new sets of variables, taking values in either the 'plus' plane or 
the 'minus' plane. The dynamics on the 'plus' plane incorporates the full dynamics 
on the (u,v) plane for A > 0, including trajectories 'at infinity' on this plane, and 
also the 'large radius' dynamics on the (u, v) plane for A < 0, again including tra- 
jectories 'at infinity' . Conversely, the dynamics on the 'minus' plane incorporates 
the full dynamics for A < and the 'large radius' dynamics for A > 0, again includ- 
ing trajectories 'at infinity'. Taken together, the dynamics on the 'plus' and 'minus' 
planes incorporates the full dynamics for either sign of A, including the trajectories 
'at infinity', which is mapped to the unit circle. The phase space for this dynamics 
is topologically a sphere. The (x + ,y + ) coordinates cover a connected open set that 
includes the whole of one 'hemisphere' while the (x_,y_) coordinates cover a con- 
nected open set that includes the whole of the other 'hemisphere'. The overlap is an 
annular region containing the unit circle on either the (x + , y + ) or y_) plane, and 
the two sets of coordinates on this region are related by 



(x+,y+) = 



^2xi + 2yi-l 



^2x{ + 2y 2 + -l 



(x+,y+) ■ 

(2.10) 
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Note that this maps the unit circle to the unit circle. Note too that the line x + = 
±l/\/2 is mapped to the line y_ = ±l/\/2, consistent with our earlier observation 
that the k — trajectories consist of the lines 2x 2 + = 1 in the 'plus' plane and 2y 2 _ = 1 
in the 'minus' plane. 



2.2 Fixed points 



Having determined the global topology, we will find it more convenient to proceed in 
terms of the new variables (x, y) defined by 



^/u 2 + v 2 + 2\A\ 



(u,v), (u,v) 



2|A| 



1 — x 2 — y 2 



[x,y) 



(2.11) 



These new variables coincide with (x +) y + ) for A > and with for A < 0, 

which means that the dynamical system on the (u, v) plane is mapped into the interior 
of the unit circle for either sign of A. We thus have two dynamical systems defined in 
the unit disc in the (x, y) plane, and we distinguish between them according to the 
sign 

r? = sgnA. (2.12) 
The equations for these two dynamical systems are 



dx X / 2 2\ ( 2 2 -i \ 

= [x -y) + {x +y - l) 

dy 1 2/2 2\ i / 2 i 2 1 \ 

- = -x (y -x) + (x +y -1) 



Y - ^ (i + v) x 2 + ( a - vP) x y 



(2.13) 



Pv - ^ (i + v) ^ + 13 (i - v) y 2 + 7^-x 2 

2 2a 



where ( is a new independent variable such that 



c = 



2 A 



1 — x 2 — y 2 



(2.14) 



The fixed points on the unit circle are clearly the same for both systems. As 
already observed, there are six such fixed points, which we consider in pairs, in each 
case giving both positions and the eigenvalues of the Jacobian matrices. Firstly, we 
have 



Two saddles 



(0,1) : {20, 



2a 



(o,-i) 



2a 



, -2(3 



(2.15) 



Secondly, there are four k = fixed points 'at infinity': 
• Two nodes: 



1 1 



TV2{ 



2a ' 2 



(2.16) 
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Figure 1: The global picture is that of a topological sphere. The equator is formed by the circle at 
infinity and the k — invariant set forms a continuous line akin to the seam of a tennis ball. Here 
we show the situation for A = 0. There are a total of 10 fixed points: six nodes (black dots) and 
four saddles (grey dots). 

• Two that are nodes for A < and saddles for A > A^: 



These two fixed points are non- hyperbolic when A = A^. 

From the global perspective, the last four fixed points are intersections of the great 
circle 'at infinity' on the spherical phase space with the continuous curve on the 
sphere defined by the combined k = trajectories of the A > and A < systems. 
If the sphere is viewed as the surface of a tennis ball then this curve can be viewed 
as its seam. This is illustrated for A = in Fig. [1] 

The structure of trajectories inside the unit circle depends on the sign of A so we 
now consider these two cases separately: 

2.2.1 A<0 

Setting r] = — 1 in (12.131) we have 




(2.17) 




(2.18) 
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• Type (i), or k = 0, fixed points at 

(x,y)=±-L-(\,\ h ) . (2.19) 
V 2\ h 

These lie within the unit circle iff A < A^. The eigenvalues of the Jacobian 
matrix at these fixed points are 

±^ {(A 2 -A 2 ), 2 (A 2 -A 2 )}. (2.20) 

This shows that these fixed points become non-hyperbolic at A = A c (where 
they coincide with the 'type (ii)' fixed points to be discussed below) and at 
A = Xh (where they coincide with fixed points 'at infinity'). 

• Type (ii) fixed points at 

(z,y) = ±— L==(2/M) . (2.21) 

v a 2 + y c 

The eigenvalues of the Jacobian matrix at these fixed points are 

± 1 {-A + - , 4 - a/A 2 - A 2 , -A - - - 4 x - a/A 2 - A 2 1 , 

2A hA /FTA 2 \ (rf-2)A V (rf-2)A V J 

(2.22) 

where 

A = 4 = . (2.23) 

x /(d-2)(10-d) 

These fixed points are saddles for A < A c , nodes for A > A > A c and foci for 
A > A. When A = A c they coincide with the type (i) fixed points. 

2.2.2 A>0 

Setting r] = 1 in (I2.13P we have 
dx 1 



(2x 2 -l) [2xy-XX h (x 2 + y 2 -l)] , 



% = ^ (2y 2 - 1) + (x 2 + y 2 -l)(P- \xy) . (2.24) 

We may again classify fixed points into two types: 
• Type (i), or k — 0, fixed points at 

(x,y) = ±-^(A,A fc ) . (2.25) 

These fixed points are inside the unit circle if A > A^. Their eigenvalues are the 
same as those of the k = fixed points for A < 0. Since A^ > A c , these fixed 
points are always hyperbolic, except when A = A^, in which case they coincide 
with fixed points 'at infinity'. 
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• Type (ii) fixed points lie on the line 2j3y = Ax, which is mapped into itself by 
the map (I2.10p that takes the exterior of the unit circle into the interior of the 
unit circle. The positions of any such fixed points are therefore given again by 
(I2.2ip . as one may verify, but these fixed points are now outside the unit circle. 

2.3 Phase portraits 

We now present a selection of global phase portraits for representative values of the 
parameters. We show the two hemispheres separately, but they fit together to form 
a topological sphere in the way described earlier. Saddle points are indicated by grey 
dots and nodes/foci by black dots. 
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Figure 2: A = 0. Hemisphere for A < 0. 



Figure 3: A = 0. Hemisphere for A > 0. 




Figure 6: X h > X > X c 
A < 0. 



Hemisphere for 



Figure 7: 
A > 0. 



Xh > X > X c . Hemisphere for 




2.3.1 Trajectories vs solutions 

For either sign of A the equations ( 12. ip are invariant under reflection through the 
origin 

(«,«)-►-(«,«). (2.26) 

It follows that this transformation takes one trajectory into another trajectory, ex- 
cept for the trajectory through the origin, which is taken into itself. In particular, 
this explains why fixed points occur in pairs, since the origin is never a fixed point 
for non-zero A. This transformation is induced by the transformation z — > —z, so 
paired trajectories correspond to domain wall solutions that are related by this dif- 
feomorphism, and z — > —z must be an isometry of the solution corresponding to the 
trajectory through the origin. 

These observations are nicely illustrated by the A = model with A < 0, for 
which all trajectories (in contrast to domain- wall solutions) are known exactly pQ. 
One finds that the two type (i) fixed points yield the metrics 

ds\ = dz 2 + e ±2/3 v^ ds 2 (e^^) . (2.27) 

These are both locally adS metrics that are foliated by (d— l)-dimensional Minkowski 
spaces. This foliation separates the entire adS spacetime into two regions^, separated 
by a Killing horizon, and each fixed point yields the metric on one of the two regions. 
The trajectory through the origin yields a solution with metric 

ds 2 = dz 2 + ^ cosh (f3 v f2\A\z S ) ds 2 (adS d -i) . (2.28) 

This is the same adS spacetime, but now foliated by (d — l)-dimensional anti-de Sitter 
spaces. Note the z —> —z isometry. Given global coordinates for the adSd-i leaves of 
the foliation, the entire adS spacetime is now covered, including both regions covered 
by the fixed-point solutions. There are also two k — 1 trajectories that yield a metric 
that is locally adS, and these correspond to two regions of adS foliated by (d — 1)- 
dimensional de-Sitter (dS) spaces. These adS examples illustrates the further point 
that a given domain wall spacetime may be represented by more than one trajectory, 
if one regards as equivalent solutions that are locally diffeomorphic. 

Another distinction between trajectories and solutions emerges when we consider 
A > 0. It is particularly instructive to consider the models with 



/ Tl 

where n is a non- negative integer; note that A = for n = and A„ < A c . For n > 1, 
solutions of the Einstein-dilaton model lift to solutions of the (d + n)-dimensional 
Einstein equations with cosmological constant A. Although the higher-dimensional 
solution will not necessarily have a domain-wall interpretation, the k = fixed point 
solutions lift to the two Minkowski-foliated locally adSd+ n metrics that are separated 



By 'adS' we mean adS and not its covering space. 
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by a Killing horizon. These metrics do not cover the Killing horizon, which becomes 
a curvature singularity of the ci-dimensional metric. However, the trajectory through 
the origin, which interpolates between the fixed points, yields a solution that lifts to 
a global covering of the adSd+ n spacetime, so this metric has a singularity at z — 0, 
so that the z — > — z transformation now interchanges two singular metrics. It is 
natural to conjecture that the same result holds for all A < A c , since the structure 
of trajectories is topologically the same as for A = 0. Because of the bifurcation at 
A = A c , there is a topology change for A > A c such that the trajectory through the 
origin now interpolates between k = fixed points 'at infinity'. 

Another instructive example is provided by the A = A/j model, for which there are 
just two k = solutions, for either sign of A, given byEI 

e XhV = \z\ exp (-\h\\A , = \z\~ x exp [~\h\\A ■ (2.30) 

Consider the A < case: then ip > for z > and if < for z < 0. These two 
possibilities correspond to the two k = trajectories interchanged by the transfor- 
mation z — > —z. Let us choose the if > trajectory with z > 0, and set A = — 2/X\ 
for simplicity. Then 

e XhV = ze?*, e XhU = z~ l e^ (A = -2/\ 2 h ) (2.31) 

and hence 

\ h a = z- 1 (z 2 - 1) . (2.32) 
We see that a — at z — 1, so that the range of the function a(z) is restricted to 

a(z) > o min = £7(1) = -— . (2.33) 

The relation of z to the affine parameter z is found by integration of 

1 ,2 

dz = ^= dz (2.34) 

'z 



We see that z ~ l^fz near 2 = 0. Thus the singularity of the metric and scalar field 
at z = is at finite affine distance. The domain wall solution in this case is singular. 



Nevertheless, it will be useful later. 



2.4 Bifurcations at infinity 

Bifurcations of families of dynamical systems occur when fixed points coincide for 
some values of the parameters, leading to a non-hyperbolic fixed point at the bi- 
furcation point in parameter space. There is a bifurcation in the A < family of 
dynamical systems when A = A c because in this case the type (i) and type (ii) fixed 
points coincide. It was shown in [1] that this is a transcritical bifurcation, in which 

5 This solution was given in pQ but only for z > 0. 
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the stability properties of two fixed points are exchanged as they pass through each 
other. 

Inspection of (12.1 7\i shows that these fixed points at infinity are non-hyperbolic 
when A = A^; this is because they then coincide with the type (i) fixed points, 
which move to infinity as A approaches A^. We now show that this leads to another 
transcritical bifurcation. We first define new variables (x, y) by 

x = ±-^= + x-y, y = ±-^= + y (2.35) 

and we define the new parameter 

s = (A - A,) . (2.36) 

We now find, for r\ = — 1, that the system is described by the equations^ 
dx 2 

— = sx + X h x 2 + —y [(1 + 2(3X h ) x-2y} + C (x, y, s) 
«<-, *h 

^ = ^y + y L (x,y) , (2.37) 

where C is a polynomial function at least cubic in its arguments, and L is a polynomial 
function that is at least linear in its arguments. 

There is a fixed point at the origin of the (x, y) plane that becomes non- hyperbolic 
at s — 0. One now considers the 'extended' system in which s becomes a variable 
subject to the trivial equation ds/d( = 0. This system has a fixed point at the 
origin and we are interested in the dynamics on the centre manifold through this 
fixed point. Because of the structure of the above equations, the equation for this 
'extended' centre manifold is simply y = through quadratic order, so the equation 
for the dynamics on the extended centre manifold to quadratic order is 

% = sx + X h x 2 . (2.38) 
dQ 

After flipping the sign of the independent variable, this equation becomes a standard 
one for a transcritical bifurcation. 



3 Flat walls 

For k = we have the dynamical system 



u = —auv + - (A — /i) (u 2 — 2A) + > 

v = —av 2 H — \uv — 2aA (3-1) 



6 The equations for x\ = 1 are slightly different but diffeomorphic and hence yield equivalent 
results. 
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The constraint is 

X 2 = e-^ [v 2 -u 2 + 2k] , (3.2) 
and it implies that we must restrict the phase space by the inequality 

t> 2 -w 2 + 2A>0. (3.3) 

Note that when v 2 — u 2 + 2 A = 0, the system reduces to the k = case of the purely 
dilaton model, so the k = trajectories of that model are the x = trajectories of 
the axion-dilaton model. 

We record here that the spacetime metric for the gauge choice (11. 6p used to obtain 
these equations is 

ds 2 = e Xr7 dz 2 + e 2 ^ds 2 (E^ d -V) . (3.4) 

3.1 Global Structure and Fixed Points 

Define the new coordinates 



The physical region v 2 — u 2 + 2A > is mapped onto the the interior of an ellipse: 

2X 2 + Y 2 < 1 . (3.6) 

This ellipse can be identified with the k = 'tennis ball seam' of the global phase 
space for constant axion trajectories described in the previous section. 

The trajectories 'at infinity' in the (u, v) plane are mapped to the intersection 
with the physical region of the orthogonal ellipse 

X 2 + 2Y 2 = 1. (3.7) 

The physical region contains (i) the region with X 2 + 2Y 2 < 1, which is the phase 
space for A > 0, and (ii) two disjoint regions with X 2 + 2Y 2 > 1, which collectively 
form the phase space for A < 0. The equations in the new variables are 

^ = \vl (1 - X 2 ) (1 - 2X 2 - Y 2 ) + i (1 - X 2 - 2Y 2 ) [\ h XY - A (l - 2X 2 )] 

^ = --nXY (1 - 2X 2 -Y 2 ) + -(l-X 2 - 2Y 2 ) [2XXY — \ h (l- Y 2 )] (3.8) 
dC. 2 2 

where ( is a new independent variable such that 



2A . . 

Observe that the boundary of the physical region, X 2 + 2Y 2 = 1 is an invariant set, 
and also that the two segments of the ellipse 2X 2 + Y 2 = 1 within the physical region 
are invariant sets. 
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3.1.1 Boundary fixed points 

There are two types of fixed point on the 2X 2 + Y 2 = 1 boundary 

• Fixed points 'at infinity': there are two fixed points on the ellipse, at (X, Y) = 
(±1,0), but these are outside the physical region. The only fixed points 'at 
infinity' that are inside the physical region are also on the boundary of the 
physical region. There are four of them, and their positions and eigenvalues are 
as follows: 

±4= (1, -1) = { — ^ , ~^ (A + A,) } (3.10) 





1 


75"' 






i 


7T- 


"7! 



V3 

Apart from the boundary fixed points 'at infinity', there is another pair of fixed 
points on the boundary provided that A ^ A/j. These lie on the intersection of 
the boundary with the line AftX = XY. They have positions and eigenvalues 

^ 1 (A, A.): i 1: A L, X2 7*~!?> (3- 11 ) 



v^TAl ' [ 2v ^TaI' v^+Al 

If A < Aft then these fixed points are on the boundary of the A > region. If 
A > A/j then there is one on each boundary of the two disjoint A < regions. 
At A = Aft these fixed points coincide with boundary fixed points 'at infinity'. 

3.1.2 Interior fixed points 

There are a pair of fixed points with positions 

(X,F) = ±A- 1 (Aft,A- /U ) . (3.12) 

where 

A = V / (A-/i) 2 + A(A- /U ) + A2. (3.13) 
These are in the allowed region 2X 2 + Y 2 < 1 provided that 

A(A-^)>A^. (3.14) 
Note that this condition requires A > /i, and that it is equivalent, for non-zero A, to 

A 2 — A 2 

/i</i c =^^. (3.15) 

These fixed points coincide with fixed points on the boundary when the inequality is 
saturated. When it is otherwise satisfied, there are two 'interior' fixed points, in the 
A > region when /i > 0, and one in each of the two A < regions when /i < 0. 
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The eigenvalues E± of the Jacobian matrix at an interior fixed point are 



E l = ±-^- Lx h + ^8fi\(X-fi) 2 -fi(8X-9i2) X\ 



E 2 = ±^- Lx h - \J 8[iX (A - fx) 2 - n (8A - 9/i) X\ 





(3.16) 



Note that 



2A 2 EiE- 



\Hi2 — — 



■fj, (A — fi) [A 2 



X\-nX] . 



(3.17) 



The last factor in this expression is necessarily positive for an interior fixed point, so 
the Poincare index is —1 for \i > and +1 for /j < 0. Index considerations therefore 
allow limit cycles only for /i < 0. However, /x < corresponds to A < 0, and in this 
case the constraint ( 13. 31) implies that ip is a monotone function, so limit cycles are 
excluded. 

3.2 Phase portraits 

We now present a selection of global phase portraits for representative values of the 
parameters. Grey dots indicate saddle points and black dots indicate nodes/foci. In 
each case the physical phase space is the interior, and boundary, of an ellipse. This 
allows for both signs of A, with A < in the upper and lower 'lobes'. 

The upper and lower segments of the elliptical boundary could be identified, in 
which case the endpoints of the two curves 'at infinity' would be joined to create a 
continuous closed curve bounding a connected A < phase space with the topology 
of a disc. This is the phase space described for /i = in [6] , where the disc was viewed 
as the 'northern' hemisphere of a 2-spher^]. However, a similar identification of the 
remaining two segments of the elliptical boundary (representing the constant axion 
trajectories for A > 0) is not possible. 



7 The analysis of [6] was for A > cosmologies but this corresponds to A < domain walls. The 
'southern' hemisphere of the 2-sphere of this reference consisted of the time-reversed trajectories for 
the same sign of A. 
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Figure 10: A < A^; \i > 0. There are no 
fixed points in the interior of the allowed re- 
gion. Saddles are grey dots, nodes black dots. 




Figure 12: A < X h ; n c < n < 0. 



Figure 11: A < X h ; \i = 0. The segments 
of the orthogonal ellipse at infinity turn into 
lines of fixed points on which there is no flow. 




Figure 13: A < A^; n < \i c < 0. There is 
now a fixed point inside the physical region. 
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Figure 14: X > \ h ; n < 0. Note the quasi- 
cyclical behaviour in the A < regions. 



Figure 15: A > A^; ^ = 0. The segments 
of the orthogonal ellipse at infinity turn into 
lines of fixed points on which there is no flow. 




3.3 Flat target space 

The fj, = case is special because the equations of motion then define a 3- dimensional 
autonomous dynamical system in terms of the variables 



0, u = &, w = x- (3.18) 



These equations are 



and the constraint is 



u = —auv + 2^ u2 ~ ^ ' (3.19) 

1, , 
w = —awv + -Xuw (3.20) 

v = -f3v 2 - —u 2 + -Xuv - —w 2 - 2(3A , (3.21) 
2a 2 2a 

v 2 - u 2 - w 2 + 2A = 4e A<T - 2 ^ • (3.22) 

For w = we recover the equations of (12. ip and the constraint (12.21) of the constant 
axion case. Moreover, fixed points of this sub-system are also fixed points of the 
3-dimensional dynamical system. 

For k = 0, we may use the constraint to eliminate the variable w from the equa- 
tions for (ii,v), which then become those of (13. ip . while the equation for w becomes 
redundant because it follows from differentiation of the k = constraint. How- 
ever, fixed points of this k = 'subsystem' are not generally fixed points of the 
3-dimensional dynamical system because there are no w ^ fixed points of the latter 
unless A = and A < (in this special case there is circle of fixed points with v = 
and u 2 + w 2 = 4ce/3\A\, but these are k ^ fixed points). What this means is that the 
fixed points of the 2-dimensional k = 'subsystem' correspond to non-trivial k = 
trajectories of the 3-dimensional system. 

It is instructive to consider the special case of A = 0. If we write 

u + iw = Re^ (3.23) 

then the equations of the 3-dimensional system reduce to the one equation ip = and 
the following two equations: 

R = -aRv , v = -av 2 + ^-R 2 - 2(3A . (3.24) 

2a 

These are just the constant axion equations (12.11) at A = with u —> R, except that 
R cannot be negative. The demonstration of section [2] that the constant axion phase 
space is globally a sphere, when allowance is made for both signs of A, now leads to 
the conclusion that the global phase space is a 3-sphere. This conclusion also holds 
when A 7^ 0, in which case ip oc sin ip, and the union of the ip = and ip = it subspaces 
may be identified with the spherical phase space of the constant axion system. 
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Figure 18: The (/z, A) plane. For parameter values inside the semi-infinite grey region A\J B there 
exist internal fixed points. The curve bounding this region from above is given by the equation 
A/i = A 2 — A?, i.e. [i — /z c . For /z < (A) these are located in the two A < regions and for /z > 
(B) they are in the A > region. The dash-dotted lines indicate the curve /zA = —2 in relation to 
the curve /z — /z c for d = 4. 



3.4 Bifurcations 

As shown by (13.171) . the pair of 'internal' fixed points are non-hyperbolic when /zA = 
A 2 — A^, which is equivalent to fi = /i c . Recall that \i < fi c was the condition for 
the 'internal' fixed points to lie within the allowed region, so when /z = \x c they have 
moved to the boundary of the allowed region, where they coincide with the boundary 
fixed points of (13. lip . Here we study the nature of the bifurcation as the curve 
H = /i c (A) is crossed in the (A,/z) plane (Fig. [TBI . 

Without loss of generality we may concentrate on the pair of fixed point in the 
Y > region, and we define new variables (X, Y) by 

X= - X Y = - Xh +X + Y. (3.25) 

^/\ 2 h + 2A 2 2A 3\\ h + 2A 2 1 ; 

The fixed point is now at the origin in the (X, Y) plane. We are interested in the 
behaviour for // ps /i c (A), so we also define 

s = . (3.26) 

The equations of the dynamical system near the origin in (X, Y, v) space take the 
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form 



together with the trivial equation s' = 0. The functions Qi and are polynomials 
in (X, Y) that are at least quadratic in these variables, with coefficients that may 
depend on s. For our purposes, the relevent terms are 

Q I = Mz*!*. + ..., ft = - W"V + - (3.28) 

The equation for the centre manifold of this system (i.e. the 'extended' centre mani- 
fold) is, to quadratic order, 



9A fc (A*-A>) (3 ' 29) 

Evidently, we must assume that A 7^ A^, since in this case there is a coincidence of 
three fixed points and the bifurcation is more complicated. The dynamics for the 
evolution on the extended centre manifold is governed by the equation 



dY _„f 2\ h (\ 2 h + 2\ 2 
d( S 3A 



sY- hK \~I > Y 2 + --- (3.30) 



From this formula it is evident that the bifurcation is transcritical. 

There is a much more complicated bifurcation at fi = 0. Again we restrict to 
A 7^ Ah,. In this case, equations (13.8j) imply that the two segments of the ellipse 
X 2 + 2Y 2 = 1 inside the physical region become lines of fixed points. The shape of 
the trajectories follows from the differential equation 

dY 2XXY — \ h (1 — Y 2 ) 



which has the solution 



dX X h XY - A (1 - 2X 2 



c{XY - \ h X f = 1 - Y 2 - 2X 2 , (3.32) 



for constant c. Trajectories inside the physical region have c > 0, and c = cor- 
responds to the boundary ellipse. All trajectories are segments of some ellipse that 
passes through the fixed points on the intersection of the boundary with the line 
XhX = XY. These fixed points are in the A < regions if A < A^ (Fig. [TT]) and in 
the A > region if A > A^ (Fig. H3|) . If A > A^, then /i c (A) > 0, so that // < \i c at the 
bifurcation. In this case each A < region is bounded by a heteroclinic cycle, and 
this leads to the quasi-cyclical A < trajectories of Fig. [TH 

As \x approaches zero, the interior fixed point approaches the line at infinity and 
merges with it at \x = 0; there this line turns into a line of fixed points. The physical 
trajectories are segments of ellipses f)3.32p emanating from the boundary fixed points 
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on either side of the x = boundary of the A > region (Fig. [T5|) . There is a 
unique ellipse, with 1/c = A 2 — A 2 , which just grazes the two boundaries at infinity 
(indicated by a dashed line in Fig. [TBI) ; it touches them precisely at the limiting 
position of the interior fixed point as p — > 0. Furthermore, the location of this fixed 
point at p = divides the lines at infinity into two regions. Concentrating on the 
Y > case, we have a line of repellers to the left of the fixed point and a line of 
attractors to the right, making this line a kind of extended saddle. As p becomes 
negative, the interior fixed point wanders into the A > region where it becomes a 
proper saddle. The drastic change in the topology of the flow is illustrated in Figs. 
[T4"t [15] and [161 If A < A^, there is no interior fixed point in the physical region as p 
goes through zero. The corresponding bifurcation is illustrated in Figs. [10] - HU At 
p = the trajectories are again given by segments of ellipses (13.321) . now emanating 
from the boundary fixed points in the A < regions. The two lines 'at infinity' are 
again lines of fixed points. There is now a universal repeller in the Y > region and 
a universal attractor in the Y < region. 



3.5 An exact axion-dilaton domain-wall solution 

The interior fixed point solution has 

(u, v) = ±A (A,, X — p), e>x = O^j « A (3.33) 

where 



Recall that the existence of an interior fixed point requires that both A — p and A/// 
be positive, and that p < p c , so that the constants A and k are real. 

Integrating the equations a = u and v — (p, we obtain the following new exact 
domain- wall solution 

a = ±X h A(z - z ) , if = ± (A - n) A(z - z ) + <^ , (3.35) 

where z and y?o are integration constants, and \ is such that 

e* IMr (x-Xo) = K, (3-36) 

where xo is a further integration constant. 

The spacetime metric for a convenient choice of (po is 

ds 2 = dp 2 + pW x ds 2 (E^- 2 )) , (3.37) 

where 

P = ±^ exp (±^A(z - z )^j . (3.38) 
This solution will play an important role in the discussion of supersymmetry to follow. 
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Although the limiting case in which /x = \x c yields a solution with constant axion, 
this case is of interest because the \x = fi c models include the d = 4 Freedman-Schwarz 
(FS) supergravity theory, which has a supersymmetric domain wall solution [T7]. In 
fact, the above solution generalizes the FS domain-wall to any d when /j = /i c . More 
significantly, it shows that the further generalization to fi < fi c requires a non-constant 
axion field. 



4 Super symmetry 

For domain wall solutions of supergravity theories it is of interest to ask of any given 
solution whether it is supersymmetric, since supersymmetry implies stability. A nec- 
essary condition for a domain wall solution to be supersymmetric is that it admit 
a Killing spinor, and this condition is also sufficient for models with only a single 
scalar field, e.g. a dilaton field, because the vanishing of the dilatino supersymmetry 
variation is then an integrability condition for the existence of a Killing spinor. An- 
other simplifying feature of single-scalar models is that both the potential V and the 
Killing spinor equation are determined by a 'superpotential' alone; in the multi-scalar 
case there is an additional dependence on the target space metric. An example of a 
single scalar supergravity model is minimal d = 3 supergravity coupled to a scalar 
multiplet. In this case the superpotential is a real function W(a) of the single scalar 
field a and the potential is given in terms of W by the formula 



V 



{W'Y-a 2 W 2 , (4.1) 



where the prime indicates a derivative with respect to a. This formula can be im- 
mediately generalized to d dimensions, as is implicit in the way it has been written. 
The analogous <i-dimensional generalization of the Killing spinor equation is 



0, (4.2) 



where is the usual covariant derivative acting on spinors. In the d = 3 super- 
gravity context this equation arises from the requirement of vanishing supersymmetry 
variation of the gravitino field, in which case e is an anticommuting spinor parameter, 
but linearity of the equation allows us to re-interpret it as a commuting spinor field, 
now on a rf-dimensional spacetime. 

The point of this generalization to arbitrary d of the notion of a supersym- 
metric domain wall solution of d = 3 supergravity is that domain walls admitting 
Killing spinors with respect to some superpotential W will be classically stable even 
when there is no underlying supergravity theory, although stability against non- 
perturbative tunnelling to some lower-energy configuration is not precluded. This 
state of affairs is described by saying that W defines a 'fake' supergravity theory. Its 
'fake' supersymmetric domain wall solutions satisfy 

/-^ = T 2ae av W , f- x a = ±2e av W , (4.3) 
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because any solution of these equations is a flat domain wall solution of the second- 
order equations for (if, a), and of the diffeomorphism constraint, such that the domain 
wall metric admits a Killing spinor. A fake supersymmetric domain-wall may be 
genuinely supersymmetric if the superpotential arises in the context of some 'genuine' 
supergravity theory with the required bosonic truncation, such that (i) the potential 
takes the form (14.11) and (ii) the condition of vanishing supersymmetry variation of 
the gravitino field reduces to (14. 2p . The conditions under which this happens for d = 5 
supergravity were investigated in [16] , and we will determine here the conditions under 
which this happens for d = 4 supergravity. When there is an underlying supergravity 
theory we shall say that the domain wall solution is 'genuinely' supersymmetric, and 
we will usually omit the adjective 'fake' that should qualify the generic case because 
this should be clear from the context. 

Only special solutions can be supersymmetric with respect to a given superpoten- 
tial W, but almost any flat domain wall determines a superpotential W with respect 
to which it is supersymmetric [H [H [5]. The proof is by construction of W. Given 
(p(z) and a choice of the function /, the first of equations H4.3[) determines W as a 
function of z. Then, given <j(z), the inverse function z(a) yields If as a function 
of cr, and one may verify that W(a) so defined yields the potential according to the 
formula (14. ip . A complex superpotential is needed for the same statement to be true 
of curved walls, but we shall restrict our attention to flat domain walls, for which 
a real superpotential suffices. The qualification 'almost' arises because the function 
z(cx) is defined only when a(z) is strictly monotonic, i.e. only when & ^ 0. This 
monotonicity condition fails for adS vacua, viewed as domain walls, because these 
have cr = 0; this is expected because adS vacua may be unstable. It also fails for 
domain walls that are asymptotic to an unstable adS vacuum because in this case 
the asymptotic solution is an accumulation point for zeros of a [14J. 

We will begin our investigation of domain wall supersymmetry by examining in 
more detail what happens at an isolated zero of &. We will confirm the 'piecewise 
supersymmetry' conclusion of [5] but the details are instructive and suggest a general 
picture that we elaborate in the context of asymptotically adS domain walls. Then 
we turn to a consideration of supersymmetry in multi-scalar models, using the axion- 
dilaton models as a 'laboratory'. The exact solution of subsection 13.51 is useful in this 
respect and what we learn from it motivates a new multi-scalar extension of the ideas 
of fake supergravity. 

4.1 Single-scalar models revisited 

The formula (14.11) implies that 

V' = AW (W" - a 2 W) . (4.4) 

It is usually concluded that stationary points of W are stationary points of V, with 
V < 0, but this assumes that W" is non-singular when W = 0. It may happen that 
a zero of W is cancelled by a pole of W" such that V is finite and non-zero. We 
shall argue here that this possibility is realized for any superpotential constructed 



28 



from a solution for which a has an isolated zero. Rather than present a general proof 
we illustrate the point with two examples, one explicit and the other implicit. The 
explicit example arises from an exact domain-wall solution of the A = A^ Einstein- 
dilaton model, for which a has a single isolated zero [1] . The implicit example arises 
in the context of asymptotically adS domain walls. 

Before proceeding we pause to discuss some consequences of the choice of the 
function / in (14. 3ft . For the choice made in (11.61) . which coincides with the choice 
made in PQ, the first-order equations ( 14. 3 j) become 

= T \ h ei x ^W , a = ±2e^ Xha W . (4.5) 

The choice / = e~ aip , which was made in [5], leads to simpler first-order equations in 
terms of an affine distance variable, which we called z in subsection 11.11 In terms of 
this affine distance variable, defined as a function of z by (II. 9p the equations (14.51) 
become 

(p = \ h W } a = -2W'. (4.6) 
where the grave accent indicates differentiation with respect to z. 



4.1.1 Branched superpotentials 

Let us first return to the domain wall solution of the A = A^ constant axion model 
with A = —2/\\ studied in subsection 12.3.11 Recall that &(z) has a zero at z — 1 for 
this solution. Our aim is to elucidate the implication of this fact for supersymmetry. 
Using the solution to construct a superpotential in the way just reviewed, one finds 
that 

W(a) = T{^P)e-^ h(T , (4.7) 
where z(a) is determined implicitly by 

e XhU = z- l e^ z ' 2 . (4.8) 
Observe that z' = \hz/(z 2 — 1), and hence that 

w '" = Tfer) e ^' (4J) 

We see that W vanishes at z — 1, i.e. at a = a m i n = cr(l). One the one hand, this is 
required by the second of equations ( 14.51) . since a — at z — 1. On the other hand, 
the absence of any stationary points of V would lead us to expect the same of W. 
This apparent contradiction is resolved by the observation that 

W" - a 2 W = T (^j^-^j e-^ Xha . (4.10) 

This diverges precisely when W' = 0, such that 



AW \W" - a 2 W] = ?-e- Xh<T = V , (4.11) 

Aft 
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a W(a) 




1 > 

&min & 

Figure 19: The branched superpotential. One branch is shown as a solid line, the other branch is 
shown as a dashed line. Both branches are defined for all a > a m . 



as expected. This illustrates the important point that stationary points of W are not 
necessarily stationary points of V. 

The superpotential of (14.71) is defined only for a > a min . This is not a problem per 
se because the solution is such that a > o~ m i n , but the minimum of a(z) is a branch 
point of the function z(o). The same is therefore true of W(a). In fact, one has 

W{a) ~ W | 0~ 0~min | 2 0~ ► 0~min 

(4.12) 

for constant Wq, implying that W(a) is double- valued with a branch point at a = 
Grain- As z varies through 1, the function W(a) switches from one branch to the 
other. In other words, the solution is 'piecewise supersymmetric'. 



4.1.2 Asymptotically adS walls 

We now move on to more general potentials V(o~). Of particular interest, and common 
occurrence in supergravity models, are potentials with stationary points of V at which 
V < 0. These yield adS vacua. Given such a stationary point at cr = 0, the potential 
has the expansion 

V =-WP + 1 2 mV + - (413 > 

where I is the adS radius and m the mass of the a-particle in this vacuum. The 
vacuum will be perturbatively stable as long as m 2 satisfies the (i-dimensional version 
of the Breitenlohner-Freedman (BF) bound [HI [19] 



m 2 >-^£=m| F . (4.14) 
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For a domain wall (flat or curved) that is asymptotic to the adS vacuum as z — > oo, 
one has [2] 

*~e-^, $~JL + ™ e -*"", (4.15) 

where v is a real positive root of the quadratic function v 2 — (d — l)v — m 2 £ 2 . The 
roots 

1 



d-l± 1) +4m 2 £ 2 



(4.16) 



are real as long as the BF bound is satisfied, and if m 2 < then both u + and v_ are 
positive. If m 2 > then only v + is positive. 

If we use the asymptotic solution (14.15p for v = v± to construct the superpotential 
with respect to which it is supersymmetric then we may assume that W > at cr = 
without loss of generality because the overall sign of W serves only to distinguish 
between walls and 'anti-walls'. One finds the superpotential 

The construction guarantees that this solves (14.11) . as one may verify, although we 
actually find the W_ superpotential this way only for m 2 > because for m 2 < 
there is no domain wall associated with W-. because v_ is not positive^. As noted 
in [1], in a somewhat different context, there is a further one-parameter family of 
perturbative solutions of (14.1 1) for which W(0) is non-zero, and a domain wall that 
is supersymmetric with respect to such a superpotential is not asymptotic to the 
adS vacuum at cr = 0. These 'other' superpotentials, Wm(&), parametrized by the 
non-zero real number M, have the expansion 



W M (a) =^M 2 + v £2 ' + aMa + -a 2 ^ M 2 + v p ' a 

For non-zero m, it is evident from the <r 3 term why M must be non-zero. For m = we 
require M^O because the M = solution then coincides with the W_ superpotential. 

To summarize, each perturbatively-stable adS vacuum is associated with two su- 
perpotentials W± that are stationary at the vacuum and a family of superpotentials 
Wm that are not. The question that we wish to address now is how the superpoten- 
tials defined at a maximum of V are related to those defined at an 'adjacent' minimum 
of V. A number of single-scalar supergravity models have a potential with a single 
adS maximum and a single adS minimum, derived from a superpotential with respect 
to which the maximum is supersymmetric and the minimum is not. The expansion of 
this superpotential about the minimum necessarily coincides with Wm for some value 



8 This generalizes an observation in [T] for the m 2 — case of constant V: in that case W- — Wq 
and W+ = Wo cosh (aa) for constant Wq but whereas both W+ and W- admit supersymmetric adS 
vacua only W+ admits a supersymmetric domain wall. 
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of M, and it appears from inspection of a few cases that it should be identified with 
W- at the maximum. In any case, any supergravity model automatically gives us at 
least one superpotential that is defined as a single-valued function for all values of 
the scalar field and hence for all values in the interval between the maximum and the 
minimum. Let us assume that we are given one such superpotential, not necessarily 
arising from a supergravity model, and ask about the other superpotentials. There 
are evidently various possibilities. It could be that each superpotential that is pertur- 
batively defined at the maximum of V continues to one that is perturbatively defined 
at the minimum. Alternatively, it could be that a pair of superpotentials defined 
perturbatively at, say, the maximum are actually two branches of a double- valued 
potential that is not defined at the minimum. It follows from recent work of Amsel 
et al. pj)] that there exist potentials for which this latter possibility is realized. We 
now present a version of their argument adapted to our needs. 

We assume (i) that V(a) has an adS maximum, at a = 0, and an adS minimum 
at a — a + but no other stationary points, (ii) that the continuation of W+, pertur- 
batively defined at the origin, leads to a function that is single valued in the interval 
[0,0+], (ii) that this function is stationary at a = a + What can we say under these 
circumstances about the continuation of the function W-. defined at the origin? Both 
u + and i/_ are real and positive at the origin, by hypothesis, so not only are both W + 
and W- positive near a = but so also are their derivatives W' + and W'_, as one sees 
from ffl~TTj) . Writing (TO) as 

y/2W' = VV + 2a 2 W 2 , (4.19) 

we see that it remains true that W + > W_ as we increase a from zero, until a zero of 
either W+ or W- is reached. Since the function V + 2a 2 W+ is strictly greater than 
the function V + 2a 2 W 2 , and since the former is zero at a = a + , there is necessarily 
a zero of W'_ at some cr_ < a + . A priori, this might be another stationary point of 
V (which would have to be some solution of W' + ' = a 2 W + since it is not a solution of 
W' + = 0, by hypothesis) but this is excluded by assumption (i), so W- is stationary 
at a point that is not a stationary point of V. As we have seen, this implies that W- 
is a double- valued function that is not defined for o > The simplest possibility 
for the second branch of this function is that it continues back to the origin, where 
it must be identified with Wm for some M. 

There is a simple interpretation of this result for the cosmological analog in which 
the potential is inverted. What was a domain wall becomes a cosmology, represented 
by the damped motion of a ball in the inverted potential —V. The adS maximum 
at the origin becomes a de-Sitter (dS) minimum and the adS minimum at o = a + 
a dS maximum. The condition that the BF bound is satisfied at the adS maximum 
becomes the condition that motion at the dS minimum is overdamped jH] , so that it 
approaches monotonically and exponentially fast with exponent u + or Suppose 
that the ball is initially at rest at the maximum and (after an infinite time) it falls 
to the minimum, exponentially fast with exponent u + . This cosmological motion 
corresponds to a domain wall interpolating between the two stationary points of the 
potential, and it determines a superpotential that (a) is single-valued in the interval 
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[0, er + ], (b) coincides with W + at the origin, and (c) is stationary at a = a + . This is the 
superpotential W + of the previous paragraph. We have assumed its existence but this 
must require a special choice of V since the ball would typically overshoot, yielding 
a superpotential with a supersymmetric adS maximum and non-supersymmetric adS 
minimum. 

Now suppose that we choose initial conditions for a ball on the other side of the 
dS minimum such that it overshoots the dS minimum at the origin with non-zero 
velocity and continues up the potential towards the maximum. It cannot 'just reach' 
the maximum (because that motion has already been accounted for) so it will either 
overshoot or it will fall back. There must be some initial conditions for which the 
latter possibility is realized. As it falls back it will either overshoot again or it will 
approach the dS minimum exponentially fast with exponent V-. Although the former 
possibility may be generic, there must exist some choice of initial conditions that 
realizes the latter. In this case the corresponding domain wall solution determines the 
double- valued potential for which the function W- defined at the origin is one branch, 
the other branch being a Wm function for some particular M (which carries the 
information about the required initial conditions in the mechanical interpretation). 
The branch point of this double-valued function is the turning point in the motion of 
the ball as it falls back towards the minimum. 

In the case that the adS maximum is unstable, the cosmological dual situation 
involves a ball that oscillates indefinitely as it falls to the minimum of the inverted 
potential. As observed in [14], this implies an accumulation of zeros of a that prevents 
a domain wall asymptotic to an unstable adS vacuum from being even piecewise 
supersymmetric. As we now see, an attempt to construct a superpotential from such 
a domain wall solution would yield a multi-valued function with an infinite number 
of branch points that accumulate at the origin, so there is no superpotential for such 
a domain wall solution that is differentiable in a neighbourhood of the adS vacuum. 

4.2 E Pluribus Unum? 

We now turn to multi-scalar models. For a generic multi-scalar model, both the po- 
tential and the Killing spinor equation must depend on both the superpotential and 
the target space metric, both of which will generically depend on all scalar fields, be- 
cause this feature is already apparent from the two-scalar models obtained by coupling 
of a chiral supermultiplet to d — 4 supergravity. Moreover, there will be additional 
requirements for supersymmetry arising from the requirement of vanishing super- 
symmetry variations for the spinor superpartners of the additional fields. However, 
a domain-wall solution of a multi-scalar model can always be viewed as a solution 
of a single-scalar model obtained by truncation after choosing 'adapted' coordinates 
on the scalar field target space for which all scalar fields but one are constant |16j . 
As the additional conditions required for supersymmetry are just the constancy of 
the 'other' scalar fields, the problem of (fake) supersymmetry for multi-scalar model 
domain-walls has been reduced to the single-scalar problem. 

However, there is a difficulty with this 'reduced' notion of fake supersymmetry (in 
addition to possible global problems arising from the local nature of the 'adapted' tar- 
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get space coordinates). It may happen, if the single-scalar truncation is inconsistent, 
that there are no nearby solutions of the single-scalar model that are also solutions 
of the multi-scalar model, in which case the fake supersymmetry of the domain wall 
as a solution of the single-scalar model is irrelevant to its stability as a solution of 
the multi-scalar model. We shall show that precisely this scenario is realized by the 
exact solution of subsection 13.51 

The relation (13.361) satisfied by the axion-domain wall solution of subsection 13.51 
suggests that we choose adapted coordinates (a, x) defined by 



e M<T = exp 



On 



J 



X ~ Xo 



exp(-^a)(x + *;)J- 1/2 , (4.20) 



where k is the constant defined in (13 .34ft . and 

1 



J = 1 + (x + k) 



Jo 



1 + -u 2 k 2 = 
4 



A(A-/i) 



\ 2 



(4.21) 



Observe that Jq > 1 for /i < /i c , and that the Jacobian of this change of variables 
is non-zero for all finite a, so the adapted coordinates are valid globally as long as 
there are no identifications on the hyperbolic target space (e.g. those arising from 
requiring periodicity in x)- I n t ne new coordinates the solution is 



x = o, 



a 



cr(z) In J 



(4.22) 



with a(z) and <p(z) as in (13.351) . This confirms that the new coordinates are indeed 
'adapted' to the solution in the sense of [16]. 

In the gauge used to find the above solution, the effective Lagrangian in the new 
variables is 

2L eff = e^-\^rh - (J/ J ) - - 2A) 



where 



and the constraint is 



A — A/ a/Jo = A? 



A 



A — jj, 



if' 2 - (J/ J ) d 2 - J- L x 2 + 2A = . 



1 ~,2 



Now observe that 

Leff = (Jo) 35 Leff 

where 



l f-i 

2 



-a<p+ iAo- 



f = e 

and that the equation of motion for /, for / as given, is 

ip 2 - 5 2 + 2A = , 



(4.23) 
(4.24) 

(4.25) 

(4.26) 
(4.27) 
(4.28) 
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which is the constraint (14.251) for \ = 0- We see that the 'adapted' truncation, 
achieved by setting \ = 0, yields a one scalar model with exponential potential. 

The above construction guarantees that the solution used in the construction of 
the 'adapted' single-scalar model is a solution of the equations of motion of this 
model. Let us verify this: the solution for (<£>,<5"), expressed in terms of A, is such 
tha|] 

, / 2A . ~ / 2A 

This is precisely the "type 1" fixed point solution of the single-scalar model described 
in pQ , which was shown in that reference to correspond to a flat domain wall solution 
that is supersymmetric with respect to the real superpotential 

W{a) = ( j 2 exp . (4.30) 

It follows that there is at least one solution of the single-scalar model that is also 
a solution of the original axion-dilaton model: namely, the solution of the latter used 
to construct the former. However, it is not guaranteed that any other solution of the 
single-scalar model will share this property; it will if the truncation is consistent, in 
the technical sense, but not necessarily if the truncation is not consistent. For the 
case in hand, the truncation is not consistent: the expansion of L e ff in powers of 
X shows that there is a linear term, which is such that the field equation for \ is 
satisfied by x = iff 

k (a 2 — 2AA//i) = . (4.31) 

This is satisfied if k — 0, but this corresponds to the limiting case in which \ = 
0, so the truncation is equivalent to setting the axion field to zero in the original 
Lagrangian; this is evidently a consistent truncation, and we implicitly used this fact 
in section [2j Otherwise, the truncation is inconsistent because it implies a condition 
on a. Using the relation between a and a at x = 0, and the constraint, we see that 
this condition states that 

a 2 = \\A\ ^ = {\-»fA\ (4.32) 

where A is the constant defined in f)3.34p . As one sees from (13.351) . this is satisfied by 
the solution we started with (as was guaranteed) but only by this solution^!. 



4.3 Axion-dilaton fake supergravity 

The above result suggests that we should extend the notion of fake supersymmetry 
to models involving at least at axion as well as a dilaton. Here we consider such an 

9 From (|4.24p it follows that A 7^ A;,, as long as fi 7^ 0, as assumed. 

10 The consistency constraint allows for the sign ofa/ip to differ from (|3.35|) but there is no solution 
corresponding to this possibility. 
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extension, using the coupling of d = 4 supergravity to a complex chiral superfield as 
a model. Consider a Lagrangian density of the form 

C =y /- det g[R- 2Gdr ■ df - V] , (4.33) 

where r is a complex scalar field taking values in a Kahler target space and 

G = d T df)C(T, f) (4.34) 

is the Kahler target space metric, expressed in terms of the Kahler potential JC. We 
take the scalar field potential V(t, f) to be given in terms of K and a holomorphic 
superpotential P(r) according to the formula 

V = l -e K [\d T P + d T K P\ 2 G- 1 - 4a 2 \P\ 2 } . (4.35) 



This generalizes to arbitrary d the standard d = 4 formula 11 !. The potential is invari- 
ant under the Kahler gauge transformation 

K^K-(f + f), P^e f P, P^e f ~P, (4.36) 

where /(r) is an arbitrary holomorphic function of r and / is its anti- holomorphic 
complex conjugate. If P = then V = 0, so only the P ^ case is of interest here, 
but for P 7^ we may choose P = 1 without loss of generality as this is a Kahler 
gauge choice. In this case we have 

V = 2 [\d T W\ 2 G- 1 - a 2 W 2 ] , (4.37) 

where 

W = e K/2 . (4.38) 
To see that the ^-dependence is correct we write 

r = x + i£(<r), (4.39) 

for some function E of the dilaton field er, and we assume that K, and hence W, is a 
function only of a, not of \. In order to get a standard kinetic term for a we must 
choose the functions £(cr) and W(a) to be related by 

W'Y /£"\ 1 

tfJ-(»J(-fJ = 2- (440) 

where a prime indicates a derivative with respect to a. Specifically, this equation 
implies that 



2Gdr -df = - 
2 



(&r) 2 + (E'r 2 (d X ) 2 . (4.41) 



The formula (14.371) now reduces to the standard d- dimensional fake- supergravity for- 
mula fl44|) . 



1 See, e.g., [5U], but note that this reference uses conventions that differ from ours. 
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For d = 4 mod 4, the Killing spinor equation may be generalized to 

e = 0, (4.42) 
where 7* is the product of the d Dirac matrices, satisfying 7^ = —1, and 



D, ~ - ^TYf K/2 [Re P + 7* Im P] r A 



A, = -~ [dfjTd T K - d^TdfK) (4.43) 

is the Kahler gauge connection; it transforms as A — > A — %d (f — f). It ensures 
invariance under Kahler gauge transformations because the spinor parameter e has 
the transformation 

e - «p (-£(/-/) 7.) e- (4-44) 
When FT is a function only of a, we have 

so it is zero for constant axion. The field strength is 

d { „A v] = -—d^ad^x, (4.46) 

but this vanishes for all domain-wall solutions, irrespective of whether the axion field 
is constant. From this we may conclude that the Kahler potential is irrelevant to 
the integrability conditions required for the existence of a Killing spinor. As these 
integrability conditions are all that we will need to consider here, we may ignore 
the Kahler gauge connection in what follows. Setting P = 1 too, we see that the 
'effective' Killing spinor equation is precisely (14.21) . This verifies that we have the 
correct d dependence in ( 14.42 p . at least for d = 4 mod 4; we will not consider here 
whether there is an appropriate generalization of (14.421) to other dimensions. 

As the potential given by (14.11) for a superpotential W that is a function only of 
a is itself only a function of a, it is manifest that the x = truncation is consistent, 
and one then recovers the fake supergravity formalism summarised at the start of 
this section. This shows how the fake supergravity formalism can be recovered from 
N=l d = 4 supergravity, and it checks the ^-dependence of our proposed extension of 
the potential and Killing spinor equations. However, for what follows it is important 
to appreciate that the formula for the potential (14.11) and the Killing spinor equation 
(14.21) are valid (at least for d = 4 mod 4) even if no truncation is made. All that 
has been assumed to derive these equations is that W depends only on a, although 
the derivation requires that the coupling of a to x then be chosen such that (14.401) 
is satisfied. If, as here, we wish to specify a particular coupling (corresponding to a 
particular choice of target space) then only certain superpotentials W will be possible. 
Specifically, our axion-dilaton model is found by choosing 

£(<t) = ±-e-2 a , (4.47) 
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in which case (14.401) reduces to 

W'Y a fW'\ 1 

w) + l{w) = r (448) 

This equation is solved b}@ 

W = W e" ltk , (4.49) 
for constant Wo, which we assume to be non-zero. We then find the potential 

V = Ae~ Xa (4.50) 

where 

A = -2/n , 2A = Wl (A 2 - A 2 ) . (4.51) 

What we have now shown is that a one-parameter subfamily of our two-parameter 
family of axion-dilaton models may be considered as 'fake axion-dilaton supergravi- 
ties', in the sense described above. This family is defined by //A = —2 and for d = 4 
there is an underlying 'genuine' supergravity theory, which includes the FS model as 
the special case for which A = 1. 

Actually, we have two separate one-parameter families: one with A > and A > A^ 
and another with A < and A < A^. In either case we have /i (A — fi) < (since 
fi = —2/ A < 0) and hence the restriction to the one-parameter family is compatible 
with the existence of the fixed point solution discussed above only if A < 0, and it 
exists in this case only if \i < fi c . In the special case that \x = /i c , which is realized by 
the FS model, we see from (I3.36j) that \ = 0, and hence that the \ = truncation 
X = is consistent. We shall exclude this case as trivial (the solution coincides with 
one of the fixed-point solutions of section [2]) so we consider only // < /i c , which is 
equivalent to A > A c when iiA = —2. To summarise, a fixed point solution with 
non-constant x exists for \x = — 2A and A < provided that 

A c < A < \ h . (4.52) 

Observe that 

for this 1-parameter family, and that W and W, of (14.301) . are identical functions 
of their arguments, as expected since the 'adapted' truncation leading to the su- 
perpotential W(a) coincides with the consistent truncation x = leading to the 
superpotential W(a) in the limit that A = A c . 

These fixed point solutions are fake supersymmetric in the sense of [5], which relies 
on the truncation in 'adapted' coordinates of [16J, but are they fake supersymmetric 
solutions of the untruncated model? We are now in a position to answer this question. 
First, we recall from [lj that the integrability conditions for the Killing spinor equation 
L2D are 

(p = ±2ae^ Xa W , a = ^2e^ Xa W' , (4.54) 



2 Othcr solutions of (|4.48[) yield a 'non-exponential' potential that is of no particular interest here. 
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and hence, for the superpotential (14.491) . 



if = ±2aW , & = ±XW . (4.55) 

Comparing with (13.351) we see that supersymmetry requires Wo = vJA) and fi = fi c 
or, equivalently, A = A c , which is outside the range (I4.52p . The fake supersymmetry 
of the limiting case for which A = A c is not surprising because x — in this case. In all 
cases for which the axion is not constant, the single-scalar truncation is inconsistent, 
the domain-wall solution is not a fake supersymmetric solution of the untruncated 
model, at least not in the sense described here. 

5 Discussion 

We have studied domain wall solutions of gravity coupled to an axion and dilaton, 
with a hyperbolic target space and exponential dilaton potential. In particular, we 
have studied the dependence on the two parameters /j ^ and A > that deter- 
mine, respectively, the target space radius and the dilaton self-coupling. Starting 
from the observation that domain wall solutions with either (i) constant dilaton or 
(ii) Minkowski 'worldvolume' are essentially determined by the trajectories of an au- 
tonomous 2-dimensional dynamical system, we have used the methods of dynamical 
systems to provide an analysis of the trajectories in these two cases, allowing for 
either sign of the potential. In case (i) our results constitute a global completion of 
previous results. In case (ii) we have both extended previous results to either sign 
potential, and provided a global completion of them. In both cases this has been 
achieved by finding appropriate new variables that allow the global picture to emerge 
naturally. We have also investigated the nature of the bifurcations that occur when 
fixed points coincide as the parameters vary. All simple bifurcations are transcritical. 
We should stress here that these dynamical systems results are equally applicable to 
domain-walls and cosmologies, but we have focused here on the domain-wall interpre- 
tation to avoid any confusion that may arise from the flip of the sign of the potential, 
and curvature k, needed to pass from one interpretation to the other. 

One new result of this paper is an exact 'scaling' solution for a flat domain wall 
with non-constant axion field. For one sign of the potential (negative for domain 
walls but positive for cosmology) this is the domain wall solution that corresponds 
to the cosmological scaling solution found in [8]. In this context, it was noticed by 
Rosseel et al. [21] that this solution is unusual in that the motion in the target space 
is not geodesic. The same observation obviously applies in the domain wall case. 

As shown in (H [TJ, [3] , all flat or adS-sliced domain wall solutions for which the 
scalar field is a strictly monotonic function are supersymmetric with respect to a 
superpotential that is determined by the solution itself. The monotonicity condition 
is violated 'maximally' by adS vacua and 'badly' by domain walls that are asymptotic 
to unstable adS vacua. Here we have investigated what happens when it is violated at 
isolated points, which is a fairly generic phenomenon. An explicit example confirmed 
the suggestion of [S] that such domain walls should be 'piecewise supersymmetric' 
with respect to a multi-valued superpotential. It also showed how this is associated 
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with stationary points of the superpotential that are not stationary points of the 
potential. A nice picture emerged from our further investigation of multi-valued 
superpotentials in which the domain wall solutions of any given potential determine 
the global branched structure of the associated superpotential, or superpotentials, 
and we used this to recover and extend results of [15] on negative potentials with 
both a fake supersymmetric maximum and a fake supersymmetric minimum. 

There is a sense in which the results on supersymmetry of domain walls in single- 
scalar models carry over, at least locally, to multi-scalar models because any domain- 
wall solution of a multi-scalar model is also a solution of the single-scalar model 
obtained by truncation in 'adapted' target-space coordinates [IB] . Obviously, one 
cannot expect to use fake supersymmetry in this sense to prove stability with respect 
to fluctuations of the truncated scalars but one might hope that it would be sufficient 
for stability against fluctuations of the one untruncated scalar. The results of this 
paper show that this hope is unfounded, in general. The problem is that the 'adapted' 
truncation is not guaranteed to be a consistent one, and if it is not consistent then 
the 'nearby' (and generically, time-dependent) solutions of the single-scalar model 
(needed to establish stability in the single-scalar context) do not lift to 'nearby' 
solutions of the original multi-scalar model. It was shown in [16J, in the context 
of d — 5 supergravity, that this problem does not arise if the original multi-scalar 
solution is supersymmetric, so only a non-supersymmetric domain-wall solutions of 
d = 5 supergravity could lead to an inconsistent adapted truncation. It would be a 
surprise if this result were not generally valid, and we conjecture that the converse is 
true too: that inconsistency of an adapted truncation is associated with a failure of 
the domain wall to be supersymmetric as a solution of an untruncated supergravity 
theory. 

We tested this conjecture on the new scaling solution mentioned above. This 
solution has the interesting feature, which we suspect is related to its non-geodesic 
character, of leading to an inconsistent one-scalar truncation. The test was made 
possible by a derivation of the d = 4 fake supergravity formalism from minimal d = 4 
supergravity coupled to a chiral supermultiplet, and an extension of this formalism 
to cover a one-parameter subfamily of our original two-parameter family of axion- 
dilaton models, which could therefore be viewed as fake axion-dilaton supergravity 
theories. The scaling solution was found to be non-supersymetric in this context, as 
the conjecture requires. 
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